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ABSTRACT 


A  detailed  numerical  model  has  been  developed  for 

characterizing  the  important  energy  transfer  processes 

operative  in  the  CO/tI2  direct  current  discharge  laser 

system.  The  model  is  based  upon  a  rate  equation  formulation 

which  includes  30  levels  of  both  CO  and  and  treats 

multi-quantum  electron -molecule  excitation  processes, 

single  quantum  vibration-vibration  exchange  and  vibration- 

translation  energy  transfer  processes  and  both  harmonic 

and  overtone  spontaneous  emission  terms.  In  this  paper 

the  time  evolution  of  the  CO  vibrational  distribution 

with  and  without  N2  is  calculated  and  the  associated  small 

signal  gain  is  predicted  and  compared  to  the  measurements 

(4) 

of  Jeffers  and  Wiswall'  .  Good  agreement  is  obtained 
between  the  predicted  and  measured  delay  times  to  maximum 
gain  in  a  pulsed  CO  laser  system  for  reasonable  assumptions 
on  the  important  system  parameters.  The  variation  in  the 
predicted  gains  and  time  delays  is  explored  for  various 
values  of  the  electron  density  and  the  translational  and 
rotational  temperature  as  a  contribution  toward  improved 
experiment  design  in  furthering  the  understanding  of  the 
phenomenology  in  the  pulsed  C0/N2  laser  system. 


Introduction 


Considerable  interest  is  being  shown  in  the  CO  laser 

(1  0) 

system  due  to  the  high  output  powers '  '  '  and  the  wide  wave- 
(3) 

length  range'  available.  As  a  result  of  this  interest, 

we  have  developed  a  numerical  modeling  code  based  upon  a 

rate  equation  approach  to  characterize  both  the  dynamic  and 

the  steady  state  (CW)  behavior  of  the  CO  system.  In  this 

paper  a  description  of  the  numerical  code  is  presented  together 

with  preliminary  results  on  modeling  the  delay  times  to 

maximum  gain  on  various  vibrational  transitions  in  a  pulsed 

CO  laser.  The  predicted  delay  times  are  compared  to  available 

(4) 

experimental  data'  '  and  the  dependence  of  the  calculated 
values  on  the  system  parameters  is  explored. 

The  dynamic  behavior  of  pulsed  CO  laser  systems  depends 
critically  on  the  various  rate  processes  which  excite  the 
molecular  vibrational  states  and  which  subsequently  produce 
rearrangement  of  this  energy  by  vibration-vibration  exchange 
processes.  These  rate  processes  are  typically  characterized 
by  short  time  scales.  On  longer  time  scales,  or  for  steady 
state  (CW)  operation,  vibration-translation  and  spontaneous 
emission  processes  become  important  in  determining  the  gain 
characteristics. 

The  molecular  laser  code  to  be  discussed  in  chis  paper 
treats  30  vibrational  levels  of  CO  and  30  levels  of  an 
additive  diatomic  gas  (N2)  with  one  inert  diluent  (He) .  All 
single  quantum  vibration-vibration  exchange  (W)  and  vibration- 
translation  (VT)  processes  are  included  for  both  CO  and  N2. 

The  particular  models  adopted  in  calculating  the  rate 
coefficients  are  described  in  the  next  section.  Multi-quantum 
electron  excitation  processes  are  included  based  upon  the 
cross  section  measurements  of  Schulz ^  and  the  calculations 
of  Chen^.  Spontaneous  emission  processes  are  included  for 
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both  the  fundamental  and  the  first  overtone  radiation  from 

(7) 

CO  as  taken  from  the  calculations  of  Young  and  Eachus 
The  values  used  in  the  code  are  given  in  the  next  section. 

In  order  to  compare  the  dynamic  behavior  of  the  CO 

vibrational  levels  with  experimental  measurements,  small 

signal  gain  is  calculated  on  each  vibration- rotation  line. 

In  the  results  section  of  this  paper  we  present  the  time 

behavior  of  the  maximum  small  signal  gain  on  each  vibrational 

transition  and  compare  these  model  calculations  with  the 

(4) 

pulsed  CO  measurements  of  Jeffers  and  Wiswall'  .  The  gain 
equations  included  in  the  code  are  also  outlined  in  the 
following  section. 

The  experimental  measurements  of  Jeffers  and  Wiswall 
were  obtained  in  a  moderately  high  pressure  -  (about  30  torr 
of  CO)  room  temperature- transverse  discharge-pul sed-laser 
system.  The  duration  of  the  pulse  discharge  was  measured 
to  be  0.6  p,sec  and  the  time  delay  until  the  peak  of  the 
maximum  gain  for  each  transition  from  6  -  5  up  to  11  -  10 
was  determined.  These  measurements  have  been  successfully 
modeled  by  the  numerical  code  presented  here  with  reasonable 
assumptions  on  the  system  parameters,  i.e.,  kinetic  and 
rotational  temperature,  electron  density,  etc.  The  variation 
in  this  nominal  calculation  to  changes  in  the  input  parameters 
has  been  studied  in  order  to  characterize  the  sensitivity  of 
the  gain  and  time  delay  predictions.  Although  the  experimental 
measurements  were  obtained  in  CO  and  CO/He  mixtures,  we  have 
also  studied  the  effect  of  added  Ng  gas  on  the  predicted  time 
delays  and  maximum  gain  curves.  These  latter  results  are 
particularly  interesting  in  that  they  point  out  the  effect  of 
anharmonic  W  processes  and  suggest  the  importance  of  gain 
measurements  with  different  concentrations  of  added  Nj  - 
keeping  the  other  system  parameters  constant  -  in  quantifying 
the  CO  pulsed  laser  phenomenology. 
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The  paper  is  concluded  with  a  discussion  on  the  agreement 
of  the  model  calculations  with  the  experimental  results. 
Suggestions  are  made  for  further  experiments  which  will  provide 
more  definitive  comparisons  with  the  current  model  parameters 
as  a  means  of  fully  understanding  the  physical  and  chemical 
processes  which  are  operative  in  the  pulsed  CO  laser  system. 
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Molecular  Laser  Model 


In  quantitatively  modeling  the  transient  and  CW  operation 
of  molecular  laser  systems,  an  understanding  must  be  obtained 
of  W,  VT,  electron-molecule  and  inert  diluent  quenching 
processes  in  addition  to  the  effects  of  emission  processes 
on  the  energy  content  of  the  vibrational  distribution. 

Although  the  current  laser  code  does  not  contain  stimulated 
emission  processes,  spontaneous  emission  terms  are  included 
as  they  can  be  a  significant  loss  mechanism  from  the  high 

/  Q) 

lying  vibrational  levels  of  CO  . 

For  a  pure  diatomic  gas  laser  system  such  as  CO,  three 
time  scales  can  be  recognized  in  the  relaxation  process. 

The  direct  excitation  of  the  CO  vibrational  levels  through 
electron-molecular  collisions  occurs  with  the  shortest  time 
scale  since  the  electron  temperature  in  direct  discharge  laser 
systems  is  high  and  the  coupling  with  vibrational  levels  is 
strong.  On  a  longer  time  scale  the  W  processes  rearrange  the 
energy  into  higher  vibrational  levels  producing  at  the  same 
time  partial  inversion  on  subsequently  higher  and  higher 

(9) 

vibrational  bands'  .  On  the  longest  time  scale,  VT  and 
radiation  processes  degrade  the  vibrational  energy  into  thermal 
energy.  It  should  be  stressed  that  these  time  scales  are  not 
entirely  separate  as  the  electron  excitation  processes  and  the 
W  processes  overlap  enough  to  have  a  significant  effect  on 
the  short  time  gain  calculated  from  the  vibrational  transitions 
(vida  infra) .  Also,  the  W  and  VT  processes  overlap,  particularly 
in  the  uppermost  levels.  In  the  presence  of  a  second  diatomic 
gas,  additional  time  scales  can  be  roughly  separated  due  to 
the  W  coupling  between  species  which  usually  occurs  on  a 
longer  time  scale  than  the  pure  gas  W  processes  but  shorter 
than  the  VT  processes  in  either  pure  gas. 
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The  total  rate  equation  used  in  the  current  molecular 
laser  model  can  be  written  as  follows: 


(1)  fiHr.  =  Y  (P  n  n-P  n  n)  +  (P '  n  n'-P'  n  n')  +  (P'  n  n'-P'  n  n') 
'  '  ,  L  sr  s  rs  r  sr  s  rs  r  sr  s  rs  r 

dt 


s=r±l 


I  |I 


)  (P  n  n  -p  n  n  ) 
sr  s  i  rs  r  m 


s=r±l  l±n  m 


*  I  1 1  _ 

s=r±l  S±R  R 


V  ,_SR  '  RS  '  . 

>  (Qd  nn-Q  n  n  ) 
sr  s  s  rs  r  «. 


+  Y  (k  nN-knN) 
,  '  sr  s  e  rs  r  e' 
s 


(Ar+l,rnr+l  “  Ar,r-lnr) 


^Ar+2,rnr+2  Ar,r-2nr^ 


in  which  nr  is  the  population  of  the  rth  vibrational  leval 
and  n  is  the  total  number  density  of  CO;  n^.  and  n'  are  the 
corresponding  densities  of  the  second  diatomic;  n"  is  the 
number  density  of  He;  p  is  the  rate  coefficient  for  the 

S  L 

VT  process 

(2)  ns  +  n  - ►  nr  +  n 

in  which  vibrational  energy  corresponding  to  E  -  E  ,  the 

til  til  ^ 

difference  between  the  energies  of  the  r  and  s  states  in 

1  m 

CO,  is  transferred  to  or  from  translational  energy;  P  is  the 

S  IT 

rate  coefficient  for  the  W  process 
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(3) 


ri  +  n.  - -  n  +  n 

si  r  m 

in  which  vibrational  energy  is  exchanged  within  CO,  the 

slight  energy  defect  being  compensated  by  translational 

SR 

energy;  and  P^r  and  Qgr  are  the  analogous  rate  coefficients 
for  the  VT  and  W  processes  involving  the  second  diatomic, 

1  •  , 

(4)  nq  +  n'  —  — »  nr  +  n* 
and 

(5)  ng  +  n*  - ^nr  +  n£  , 

k  is  the  rate  coefficient  for  excitation  of  CO  from  the  s 
sr  th 

to  the  r  by  electrons  of  number  density  N  ;  and  the  A 

w 

coefficients  are  spontaneous  emission  rate  coefficients. 

In  «sn  earlier  work^10,11^,  we  concluded  that  in  the  pure 

gas  case  only  single  quantum  vibrational  transfer  processes 

need  to  be  considered.  However,  for  mixtures  of  gases  with 

greatly  differing  vibrational  spacing  there  is  a  possibility 

that  accidental  resonant  vibration-vibration  exchange  might 

involve  a  double  quantum  transition  in  one  of  the  species. 

In  general,  for  the  lower  vibrational  levels,  these  processes 

will  be  slow  by  comparison  to  single  quantum  processes  due 

to  the  form  of  the  rate  expressions^^' and  we  will  restrict 

ourselves  in  this  report  to  systems  wherein  only  single  quantum 

processes  are  important.  It  is  important  to  point  out,  however, 

that  for  W  processes  involving  higher  lying  vibrational 

transitions,  multiquantum  processes  can  have  larger  rate 

(13) 

coefficients  than  single  quantum  processes'  '  .  Due  to  the 
low  density  of  the  uppermost  vibrational  levels,  the  rates 
associated  with  these  processes  are  usually  small  by  comparison 
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to  processes  involving  exchange  with  the  low  vibrational 
levels  and  thus,  the  approximation  of  single  quantum  W 
exchange  processes  remains  good  for  most  practical  model 
applications . 

The  full  equation  as  shown  in  equation  (1)  can  be 
somewhat  simplified  through  the  use  of  the  detailed  balance 
relationships,  which  may  be  written 


(6) 


*  .  ,  =  pr  +  l  r  e“*Er  +  1  "  Er^kT. 
rf  r  +  1  r  +  1,  r  * 


Qr,  r  +  1  "  Qr  +  1,  re 


“(Er  +  1  _  Er)/kT, 


ml  =  plm  e  -  <E,  +  E,  -  Er  -  EJ/kT, 
rs  sr 


RS  SR  -(Es  +  Fs  -  Er  -  F  )/kT, 

=  Q _ e 


‘rs  sr 


and 


k 


rs 


*  k 


sr 


-  e  )  At 
r  , 


where  Fj  represents  the  vibrational  energy  of  the  i-th  level 
of  N2  above  the  zeroth  level. 

The  equations  as  shown  in  (1)  above  including  the  detailed 
balance  relationships  shown  in  (6)  and  using  the  rate  coefficients 
outlined  in  the  following  subsections  have  been  numerically 
integrated  by  the  Runge-Kutta-Merson  method  for  a  set  of  30  levels 
of  CO  and  30  levels  of  N2.  This  numerical  technique  has  been 
used  extensively  for  the  integration  of  large  sets  of  ccupled 
non-linear  ordinary  differential  equations  such  as  are  found 
when  modeling  coupled  chemical  reactions.  A  particular  feature 
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(14) 

of  the  code  as  taken  from  the  work  of  Keneshea'  '  is  the  use 
of  a  version  of  the  steady  state  (S-S)  approximation  for  those 
equations  which  satisfy  an  appropriate  criterion.  In  a 
coupled  set  of  rate  equations  such  as  occur  for  vibrational 
relaxation,  the  equation  can  be  written 


(7) 


Epr  "ERr  ‘  nr 


where  the  F  and  terms  may  be  complicated  functions 

of  all  the  other  concentrations  and  rate  coefficients.  Under 
conditions  where  the  magnitudes  of  both  £Fr  and  the  term 
^Rr  .  nr  are  large  but  their  difference  very  small,  the 
time  progress  in  the  normal  Runge-Kutta  integration  is  very 
slow,  i.e.  the  set  of  equations  become  "stiff".  In  the  Keneshea 
version  of  the  code  as  used  here,  each  equation  is  checked  to 
determine  if  it  satisfies  a  S-S  criterion,  meaning  the 
magnitude  of  the  £  Fr  term  is  large  but  the  derivative  dnr/dt 
is  small.  When  this  criterion  is  satisfied  the  equation  is 
iterated  over  the  At  time  step  using  an  algebraic  algorithm  as 
shown  below. 


For  conditions  under  which  F  and  Rr  *  nr  are  koth 
large,  rather  than  use  the  canonical  S-S  assumption 


(8) 


0 


which  leads  to 


(9) 


t 


the  more  realistic  assumption  is  made  that  the  £ Fr  and  S  Rr 
terms  are  effectively  constant  over  the  At  step  leading  to  a 
first  order  ordinary  differential  equation.  The  solution  is 
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•  At 

ZRr  ' 


This  solution  is  used  for  all  equations  satisfying  the  steady 
state  criterion  and  is  iterated  over  a  it  step  together  with 
the  other  equations  which  continue  to  be  integrated  by  the 
Runge-Kutta-Merson  technique.  When  an  equation  ceases  to 
satisy  S-S,  it  is  automatically  returned  to  the  R-K-M 
subroutine  for  integration. 

By  the  use  of  this  S-S  algorithm,  the  time  step  in  the 
integration  can  be  kept  large  while  not  sacrificing  either 
stability  or  accuracy  in  the  coupled  set  of  equations. 


(10) 


nr  ( t  +  At)  = 


=  (nr<t>  -  -jhrr ) 
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Vibration-Vibration  Exchange  Rate  Coefficients 


In  deriving  a  set  of  W  rate  coefficients  for  use  in 

the  exchange  processes  involving  pure  CO  collisions,  use 

was  made  of  the  recent  data  of  Hancock  and  Smith^3'^*^. 

This  data  has  been  interpreted  as  involving  long  range  forces 

for  near  resonant  W  transfer  while  the  non-resonant  W 

exchange  processes  are  dominated  by  short  range  forces.  In 

our  previous  calculations  on  W  rate  processes we 

have  only  included  the  effects  of  short  range  forces  through 

(12) 

a  modified  Rapp-Englander  Golden  W  model'  '  .  The  expression 
used  for  these  rate  coefficients  is 


11,  'S  -  508-iy/V0/ V^/3.5/6^ 

* { exp  -1.32(LEd)2/3  £  1/3 J 

where 

p,  =  reduced  mass  in  molar  units 
L  =  exponential  potential  parameter  (A) 

o 

n  =  hard  sphere  diameter  (A) 

T  =  kinetic  temperature  (°K) 

Ed  =  energy  defect  between  vibrational  quanta  being 
exchanged  (cm-1) 

and 


U  =  matrix  element  for  the  r 


n  vibrational  transition 


The  matrix  elements  used  are  modified  Harmonic  Oscillator  (HO) 
(17) 

matrix  elements '  '  in  which  the  actual  energy  spacing  between 

vibrational  units  is  used  instead  of  the  constant  spacing  assumed 
in  the  H-0  model .  Very  close  agreement  has  been  found  for  these 
matrix  elements  by  comparison  to  Morse  Oscillator  results . 
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The  vibrational  level  dependence  of  these  rate  coefficients 
does  not  match  the  dependence  as  determined  by  Hancock  and 
Smith  due  to  the  neglect  of  long  range  forces  in  the  W  model 
as  shown  in  equation  (11) .  We  have  therefore  developed  an 
analytical  expression  which  when  applied  to  the  above  equation 
produces  a  forced  fit  to  the  experimental  data  at  300°K.  The 
resulting  rate  coefficients  for  CO  are  shown  in  Figure  1  in 
comparison  to  the  available  experimental  data.  The  temperature 
dependence  of  the  rate  coefficients  is  assumed  to  be  given  by 
equation  (11),  i.e.,  by  the  short  range  force  expression. 

For  the  CO-N2  VV  exchange  rate  coefficients  in  which  the 
anharmonicity  is  considerably  greater  than  in  the  pure  gas, 
the  short  range  force  expression  was  directly  used.  This  is 
expected  to  be  a  reasonable  approximation  for  this  collision 
case. 
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Figure  1: 


Vibration  -  Vibration  Exchange  Rate  Coefficients  for  the 
reaction  CO(v)  +  CO ( J)  CO (v-1 )  +  CO(J+l)  at  300  °K. 


Vibration-Translation  Rate  Coefficients 


Vibration-translation  rate  coefficients  have  been  taken 

(18) 

from  measured  relaxation  times'  ’  using  a  harmonic  oscillator 
(19) 

model'  .  The  expression  used  is 


(12) 


where  pr  is  the  measured  relaxation  time.  Scaling  with 
vibrational  level  is  given  by  the  square  of  the  vibrational 
matrix  elements. 

Electron-Molecule  Rate  Coefficients 

The  electron  rate  coefficients  are  calculated  from 
cross  section  data  of  Schulz  ^  and  Chen^  .  The  rate 
coefficients  are  defined  in  terms  of  the  cross  section, 

Q(V) ,  as 

(13)  k(T)  =  J  Q(v)  v  f(v)  dv 

where  f(v)  is  the  Maxwell-Boltzmann  distribution  of  velocities. 
By  transforming  f(v)dv  to  a  distribution  of  speeds,  equation  (13) 
can  be  written 


3/2  . 

(14)  K(T)  =  4it  J  0(v)v»  exp 

Equation  (14)  can  be  easily  changed  to  the  energy  representation, 
yielding  the  following  rate  coefficient: 


-mv* 


2kT 


dv, 


13 


3/ 2 

v2  r 

ns,  km  -  ®tt  (  2Jt  )  ( 

— J  Q(E)  E  exp 

”E/kT 

2 

where  E  =  %mv  ,  m  being  the  reduced  mass  of  the  collision 
(essentially  the  electron  mass) . 

The  available  data  on  electron  excitation  cross  sections 
for  Nj  and  can  be  fit  over  a  limited  electron  energy 

range  E^—  E^,  where  >  E1  by 

(16)  Q(E)  «  AE  +  B 

where  A  and  B  are  constants  determined  by  the  end  point 
values  E^,  Q1  and  E2  ,  Q2. 

Over  the  energy  range  E^ — E2  =  AE12,  equation  (15)  can 

be  integrated,  yielding 


5i  "Vi  ' 

•v. 

~E/kT .. 

(17)  k(T)  =  8n(2nkT) 

(nt) 

e  (kT) 

Ea 

• 

Ea* 

+AEa 

and 

lI3 

m 

(2kTA+B)  (E+kT) 


E, 


(18) 


k(T)  =7  k(T) 

t-i 


Equations  (17)  and  (18)  have  been  solved  numerically  to 
yield  excitation  rate  coefficients  as  a  function  of  the  electron 
temperature  using  input  cross  section  data. 
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The  rate  coefficients  for  the  excitation  of  vibrational 

levels  of  N2  and  CO  by  electrons  in  the  electron  temperature 

range  1000 °K  to  30,000°K  using  the  measured  cross  section 
(5) 

data  of  Schulz'  are  shown  in  Figures  2  and  3,  respectively. 

The  complete  numerical  results  for  these  cases  are  available 
elsewhere ,  together  with  the  cross  section  data  used  in 
the  calculation,  in  general,  the  error  limits  on  the  calculation 
rate  coefficients  are  about  ±  20%. 

The  spacing  of  the  electron  excitation  rate  coefficients 
for  progressively  larger  quantum  changes  reflects  the  negative 
ion  mechanism  responsible  for  the  excitation  process.  In 
both  N,  and  CO,  the  negative  ions  are  unstable  by  1.89  eV^'2^ 
and  1.3  eVv  ,  respectively,  corresponding  to  about  the  6th 
vibrational  level  of  each  molecule.  The  spacing  of  the  rate 
coefficients  indicates  a  rather  close  coupling  between  levels 
near  v  =  6  while  quantum  changes  greater  than  £v  =  6  show  a 
considerably  reduced  coupling,  i.e.  a  wider  spacing. 

The  excitation  rate  coefficients  calculated  from  the 
theoretical  cross  sections  of  Chen^  are  tabulated  together 
with  the  corresponding  cross  sections  elsewhere ^20^ .  The 
general  agreement  between  the  calculated  and  measured  values 
is  within  the  accuracy  of  the  overall  rate  coefficients, 
i.e.  20%. 

By  observing  the  calculated  rate  coefficients  derived 
from  Chen's  cross  sections,  some  indication  of  the  scaling 
necessary  with  increasing  vibrational  level  can  be  obtained. 
Figures  4  and  5  show  this  relationship  for  &v  =  1  and  &v  =  5 
processes,  respectively.  No  trend  in  the  calculated  values 
are  noted  with  increasing  vibrational  level  in  the  electron 
temperature  range  between  10,000°K  and  30,000°K,  therefore, 
to  within  the  accuracy  of  these  calculations  all  excitation 
processes  with  the  same  &v  change  can  be  assigned  the  same 
value  as  the  v  =  0  excitation  case  at  any  given  electron 

temperature. 
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EXCITATION  RATE  COEFFICIENT  (cc/ molecult -sic) 


EXCITATION  RATE  COEFFICIENT  (cc/molecul«-sec) 


EXCITATION  RATE  COEFFICIENT  (cc/molecule-**$) 


Figure  4: 

VIBRATIONAL  EXCITATION  RATE  COEFFICIENTS 
FOR  N. (v-^v1)  BY  ELECTRONS 


ELECTRON  TEMPERATURE  ^Kxltf3) 


It  is  also  very  important  to  note  that  above  about 

10,000°K  the  electron  excitation  rate  coefficients  are 

effectively  independent  of  temperature.  This  fact  means 

that  in  calculating  the  vibrational  distributions  in  a  discharge 

environment  in  which  the  electron  temperature  is  greater  than 

this  value,  no  account  has  to  be  taken  of  the  non- Boltzmann 

122) 

electron  distributions  calculated  to  be  present'  .  This 
is  not  the  case  for  excitation  of  electronic  states  since 
these  processes  will  be  very  sensitive  to  the  high  energy 
electron  tail. 


Spontaneous  Emission  Rate  Coefficients 


As  part  of  the  modeling  of  molecular  lasers,  the 
spontaneous  emission  from  vibrational  transitions  in  the 
infrared  must  be  characterized.  This  serves  two  purposes; 
firstly,  the  emission  from  higher  lying  vibrational  levels 
can  be  a  significant  energy  loss  mechanism  in  distorting  the 
vibrational  distribution  and  secondly,  in  evaluating  the 
model  laser  code  through  side-light  emission  the  code  must 
be  able  to  generate  artificial  spectra. 

Young  and  Eachus^ have  evaluated  4  cubic  dipole  moment 
functions  and  have  concluded  that 

(19)  M(r)  =  ±0.112  +  3.11(r  -rfi)  -  0.15(r  -  rg) 2-2 . 36 (r-re) 3 
gave  the  best  fit  to  measured  relative  rotational  line 

(7) 

intensities  in  CO.  Using  the  integrated  band  intensities'  ' 
for  the  fundamental  and  overtone  bands  of 

“fundamental  "  261  amagat-1™-2 

and 

“overtone  =  1*88  ama*a t ^  cm"2  ' 

the  following  table  of  Einstein  Coefficients  has  been 
generated  and  incorporated  into  the  laser  model  code.  Current 
work  is  concentrated  on  a  subroutine  to  generate  the  entire 
artificial  spectrum  in  the  vibration-rotation  bands  including 
an  appropriate  slit  function. 


21 


Table  1. 


Relative  Einstein  Coefficients  for  Spontaneous  Emission  for  CO 


upper  vibrational 


level 


fundamental 


overtone 


1 


1.00 


2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 
29 


1.93 
2.78 

3.53 
4.21 

4.81 

5.35 
5.74 
6.23 
6.58 
6.88 
7.14 

7.35 
7.60 

7.73 
7.83 
7.90 

7.93 

7.94 
7.83 
7.69 

7.53 

7.35 
7.16 
7.00 

6.74 
6.51 
6.27 

5.82 


1.00 

3.00 

5.96 

9.89 

14.74 

20.48 

27.03 

34.36 

42.43 

51.17 

60.54 

70.79 

81.62 

92.99 

104.82 

117.05 

131.84 

147.24 

160.27 

173.80 

187.53 

201.19 

214.94 

228.60 

242.13 

255.45 

268.52 

281.40 
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Small  Signal  Gain 


In  order  to  model  the  behavior  of  pulsed  and  CW  laser 

operation  we  have  developed  a  numerical  routine  to  calculate 

the  gain  associated  with  a  particular  vibration-rotation 

line.  The  equation  for  small  signal  gain  can  be  found  in 

(23  24) 

several  references  '  '  to  be 


where  Nv,  and  Nv  are  the  vibrational  populations  of  the 
lower  and  upper  levels  respectively,  Sj  =  J  +  1  for  p  branch 

transitions,  S_  =  J  for  R  branch  transitions  (J  is  the 

U 

rotational  quantum  number  of  the  upper  level) ,  F(J)  is  the 

energy  of  the  Jth  rotational  level  and  is  approximated  by 

2 

BJ(J  +  1)  ,  and  R^,  is  the  vibrational  matrix  element  for 
the  transition  v  to  V . 

Since  the  current  application  for  the  pulsed  laser  code 
is  to  systems  characterized  by  rotational  temperatures,  we 
can  either  calculate  the  gain  associated  with  any  particular 
laser  line,  i.e.  select  both  a  v  and  J,  or  we  can  calculate 
the  J  value  for  maximum  gain  and  thus,  generate  the  maximum 
gain  associated  with  a  given  vibrational  transition.  The 
results  as  shown  in  the  next  section  represent  the  maximum 
gain  on  each  vibrational  band. 
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Results 


The  laser  code  outlined  above  has  been  applied  to  the 

(4\ 

pulsed  CO  laser  experiments  of  Jeffers  and  Wiswall'  . 

These  experiments  were  performed  at  roughly  30  torr  on  room 

temperature  CO  with  electron  pulse  times  of  about  0.6  ^sec. 

The  laser  code  was  thus  applied  to  an  ambient  CO  vibrational 

distribution  assuming  various  square  wave  pulse  shape  electron 

concentrations.  Model  calculations  were  made  with  time 

dependent  kinetic  and  rotational  temperature  in  which  the 

time  dependent  kinetic  temperature  produces  time  dependent 

rate  coefficients  for  the  W  and  VT  processes  while  the  time 

varying  rotational  temperature  enters  directly  into  the  gain 

calculations.  The  closest  agreement  with  the  experimental 

results  was  obtained  assuming  that  the  electron  density  was 
13 

10  electrons/cc,  the  kinetic  and  rotational  temperatures 
being  constant  at  T  =  300 °K,  and  with  an  electron  pulse  width 
of  0.6  p,sec.  The  results  of  these  calculations  are  displayed 
in  Figures  6  through  9. 

Figure  6  shows  the  time  evolution  of  the  initial  ambient 

13 

vibrational  distribution  of:  CO  due  to  10  electrons/cc  with 
a  Maxwellian  distribution  at  20,000°K.  As  indicated  in  the 
previous  section  on  electron  excitation  rate  coefficients, 
the  electrons  tend  to  preferentially  excite  vibrational  levels 
from  v  =  4  to  v  =  8  due  to  the  negative  ion  excitation 
mechanism.  This  behavior  is  clearly  indicated  in  the  figure. 
At  longer  times,  the  preferential  excitation  is  smoothed  out 
by  excitation  to  higher  lying  levels  and  by  W  processes 
which  become  significant  for  times  of  the  order  of  the  pulse 
duration.  The  affect  on  the  vibrational  distribution  due  to 
the  electron  pulse  of  W  processes,  electron  concentration, 
and  pulse  width  is  given  in  Figure  7. 
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Figure  6:  Evolution  of  the  CO  vibrational  distribution  due 
to  hot  electrons  in  a  pulsed  CO  laser. 


uutteu  curve:  witnout  W  processes 
dashed  curve:  T  function  time 
upper  solid  curve:  nominal  case 
lower  solid  curve:  Ng  =  3.012e/cc 
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Figure  8:  Evolution  of  the  CO  vibrational  distribution  due 
to  W  processen  following  a  0.6  jisec  pulse  for 
a  pulsed  CO  laser. 
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Figure  9:  Time  sequence  of  maximum  gain  for  the  fundamental 

transitions  in  a  pulsed  CO  laser.  Arrows  represent 
data  of  Jeffers  and  Wiswall. 
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Without  W  processes,  as  shown  in  Figure  7,  the 
vibrational  distribution  is  smoother,  i.e.  less  structured, 
than  with  W  processes  due  to  the  redistribution  of  energy 
from  low  vibrational  levels  to  higher  levels  as  a  result  of 
the  anharmonicity  in  the  vibrational  spacing.  As  the  kinetic 
temperature  is  changed  initially  from  500°K  to  300 °K  with  a 
time  constant  of  2  ^sec,  the  redistribution  of  vibrational 
energy  is  further  enhanced  and  the  vibrational  distribution 
becomes  nearly  inverted  between  levels  v  =  5  and  v  s=  7. 

When  the  electron  density  is  reduced  the  energy  supplied 
into  the  CO  distribution  is  decreased  as  is  clearly  shown  in 
the  figure.  The  effect  of  these  variations  on  the  subsequent 
maximum  gain  predictions  is  shown  in  the  later  figures  of 
this  section. 

Following  the  electron  pulse,  the  vibrational  distribution 
is  redistributed  entirely  due  to  W  processes.  The  evolution 
in  the  vibrational  distribution  of  CO  is  thus  determined  as 
shown  in  Figure  8.  Clearly,  energy  flows  up  the  vibrational 
structure  as  expected  from  anharmonicity  arguments,  however, 
due  to  the  preferential  excitation  by  electrons, vibrational 
energy  also  flows  into  the  v  =  1  level  from  levels  between 
about  v  =  4  to  v  =  8.  This  "reverse*  flow  of  vibrational 
energy  tends  to  reduce  the  possibility  of  obtaining  laser 
action  in  the  low  vibrational  levels  on  time  scales  greater 
than  the  pulse  length  without  the  presence  of  an  additive  gas 
which  would  preferentially  reduce  the  excitation  in  these  low 
levels.  This  point  is  discussed  further  in  the  discussion 
section. 

As  shown  in  Figure  8,  inversion  is  predicted  to  move  up 
the  vibrational  structure  creating  a  favorable  condition  for 
laser  gain  on  higher  vibrational  transitions  with  increasing 
time,  as  has  been  experimentally  observed^' .  The  actual 
maximum  gain  predicted  from  this  "nominal"  case  is  presented 
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in  Figure  9  where  the  measured  time  delays  of  Jeffers  and 

Wiswall  are  indicated  by  arrows.  The  agreement  with  the  time 

delays  to  maximum  gain  for  the  7-6,  9-8,  10-9  and  11-10 

transitions  are  very  good  but  the  6-5  transition  is  predicted 

to  be  about  1  p,sec  sooner  than  observed  and  the  8-7  transition 

is  predicted  to  be  about  1  p,sec  later  than  observed.  Both  of 

these  transitions  are  very  sensitive  to  the  electron  excitation  rate 

coefficients  since  they  are  predicted  to  reach  maximum  gain 

only  a  few  p,sec  after  the  end  of  the  electron  pulse.  We 

have  not  yet  tried  to  vary  the  electron  rate  coefficients  to 

improve  the  short  time  scale  predictions  since  the  agreement 

at  the  longer  time  scales  is  good.  In  the  following  figures 

the  variation  in  these  "nominal"  gain  predictions  are  studied 

for  several  different  values  of  the  important  system  variables. 

Each  subsequent  figure  has  as  reference  points  the  delay  time 
and  maximum  gain  associated  with  the  "nominal"  predictions 
together  with  arrows  indicating  the  extent  of  variation  from 
these  "nominal"  predictions. 

Figure  10  shows  the  greatly  decreased  gains  predicted 

12 

when  the  electron  density  is  reduced  to  3  x  10  electrons/cc 
while  keeping  the  other  system  parameters  constant.  The  results 
for  this  calculation  are  as  expected;  since  the  electron  density 
is  lower  the  energy  input  into  the  CO  vibrational  levels  is 
lower  and  the  maximum  gain  is  lower  (see  equation  20)  .  In 
addition,  since  the  population  of  the  vibrational  levels  is 
lower,  the  rate  at  which  energy  is  redistributed  due  to  W 
processes  is  slower  than  the  "nominal"  case  and  the  delay  times 
are  thus  longer.  Note  the  significant  effect  of  a  lower 
electron  density  on  the  maximum  gain  predicted  from  the  upper 
vibrational  transitions. 

Since  Jeffers  and  Wiswall  indicated  that  the  kinetic  and 
rotational  temperature  in  their  pulsed  experiments  were  increased 
due  to  energy  supplied  by  the  electron  pulse,  we  altered  the 
code  to  permit  use  of  a  time  varying  kinetic  and  rotational 
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Figure  10:  Time  sequence  of  maximum  gain  for  the  fundamental 
transitions  in  a  pulsed  CO  laser.  The  points 
refer  to  "nominal"  pulsed  CO  laser  predictions. 


temperature  having  a  time  constant  of  2  ^sec.  This  time 
constant  appears  to  be  consistant  with  the  experimental 
observations  although  quantitative  results  on  this  effect 
still  need  to  be  obtained.  We  have  arbitrarily  assumed  that 
the  temperature  function  goes  from  500 °K  to  300 °K  but  there 
is  no  concrete  evidence  to  support  this  particular  assumption. 
Two  features  of  the  calculation  as  displayed  in  Figure  11 
are  important  to  note;  first  that  the  time  delays  are  uniformly 
shortened  and  that  the  maximum  gain  is  decreased  over  the 
"nominal"  case.  The  fact  that  the  time  delays  are  shortened 
is  a  result  of  the  increased  W  rate  coefficients  at  the 
higher  kinetic  temperatures  while  the  decreased  gains  are  a 
result  of  the  increased  temperature  on  both  the  gain  calculation 
and  the  W  processes.  At  higher  kinetic  temperatures,  the 
exothermicity  involved  in  transferring  the  larger  (and  lower) 
vibrational  quanta  into  the  smaller  (and  higher)  vibrational 
levels  is  reduced .  This  produces  less  of  a  tendency 
toward  vibrational  inversion  and  reduces  the  expected  gain. 

In  addition,  since  the  rotational  temperature  is  also  higher, 
the  gain  is  directly  reduced  by  spreading  out  the  vibrational 
population  over  more  J  states  with  a  subsequent  reduction  in 
gain. 


To  test  the  sensitivity  of  the  predicted  gain  characteristics 
and  time  delays  to  the  electron  pulse  width,  we  increased  the 
“nominal"  pulse  width  from  0.6  ^sec  to  0.8  (j,sec  while  keeping 
the  other  parameters  constant.  The  resulting  maximum  gain 
curves  are  shown  in  Figure  12.  The  time  delays  are  uniformly 
shortened  and  the  gains  for  low  levels  are  decreased  while 
for  higher  transitions  the  gains  are  increased.  The  decreased 
time  delays  are  a  result  of  the  increased  vibrational  populations 
and  subsequent  increased  W  rates.  The  gain  results  for  low 
levels  are  a  result  of  producing  a  smaller  inversion  due  to  the 
increased  effect  of  electrons  (being  present  for  0.8  y,sec  as 
opposed  to  0.6  p,sec)  .  This  point  is  made  clearer  by  reference 
to  Figure  4  where  the  vibrational  distribution  due  to  electrons 
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MAXIMUM  GAIN  (CM1) 


Figure  11:  Time  sequence  of  maximum  gain  for  the  fundamental 
transitions  in  a  pulsed  CO  laser  with  variable 
kinetic  and  rotational  temperature.  The  points 
refer  to  "nominal"  pulsed  CO  laser  predictions. 
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Figure  12:  Time  sequence  of  maximum  gain  for  the  fundamental 


transitions  in  a  pulsed  CO  laser.  The  points 
refer  to  "nominal"  pulsed  CO  laser  predictions. 


both  with  and  without  W  processes  is  shown  for  a  0.6  ^sec 
pulse.  From  this  figure  it  can  be  noted  that  the  electrons 
tend  to  give  a  smoother  vibrational  distribution  without 
than  with  W  processes.  Therefore,  the  distribution  due  to 
0.8  p,sec  electrons  will  tend  to  be  smoother  on  time  scales 
of  a  few  (j,sec  than  for  only  0.6  jj,sec  pulse  widths  even  though 
the  total  energy  in  the  distribution  is  greater.  This  latter 
fact  is  reflected  in  the  longer  time  scales  of  figure  12 
(and  higher  transitions)  where  the  gain  is  seen  to  increase 
over  the  "nominal"  case.  Thus,  after  a  time  sufficient  for 
the  initial  distribution  to  be  destroyed,  the.  W  processes 
rearrange  the  energy  consistant  with  the  ratio  of  the 
anharmonicity  to  kT. 

By  contrast  to  this  case,  we  calculated  the  effect  of 
added  N2  on  the  pulsed  CO  laser  performance.  In  a  broad 
sense,  since  N2  is  an  additional  sink  for  electron  energy 
and  since  N2  will  tend  to  preferentially  pump  CO  by  vibration 
vibration  exchange  processes  due  to  the  smaller  spacing  of 
CO,  we  might  expect  to  find  increased  laser  gains  by  the 
addition  of  to  the  system.  This  expectation  is  not 
observed  as  is  shown  in  figures  13  through  15. 


The  predicted  maximum  gains  are  shown  in  figure  13  for 
the  case  in  which  the  N2  concentration  is  twice  the  CO 
concentration.  We  note  that  the  delay  times  are  shortened 
especially  for  the  higher  order  transitions  which  the  gains 
are  uniformly  decreased.  The  reasons  for  this  behavior  are 
a  direct  result  of  the  anharmonic  character  of  W  processes. 
Although  N2  picks  up  a  considerable  amount  of  energy  from  the 
electrons  during  the  pulse  time,  the  important  point  is  the 
mechanism  by  which  N2  vibrational  energy  couples  to  the  CO 
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Figure  13:  Time  sequence  of  maximum  gain  for  the  fundamental 
transitions  in  a  pulsed  C0/N2  laser.  The  points 
refer  to  "nominal"  pulsed  CO  laser  predictions. 


vibrational  population,  i.e.  how  does  vibrational  energy  flow 
between  Nj  and  CO.  Due  to  the  layer  vibrational  spacing  of 
Nj  and  the  fact  that  the  W  rate  coefficients  generally 
increase  the  closer  the  resonance  in  the  vibrational  quanta 

being  exchanged,  vibrational  energy  from  N2  tends  to  flow  into 
CO  via  the  low  levels  where  the  resonance  is  closest.  Thus, 
although  CO  picks  up  a  large  amount  of  energy  from  N2  it  does 
this  in  a  way  which  decreases  the  inversion  between  levels 
and  therefore  produces  decreased  gains.  However,  since  the 
low  levels  of  CO  have  a  layer  population  than  the  "nominal" 
case,  W  rates  in  CO  are  higher  and  hence  the  time  delays  are 
shorter. 

A  further  understanding  of  this  N2-CO  transfer  can  be 

seen  by  reference  to  figures  14  and  15  where  the  vibrational 

distribution  due  to  the  electrons  and  subsequent  W  relaxation 

are  displayed.  Since  the  W  transfer  from  N~  to  CO  is 

-1  ^ 

accompanied  by  several  hundred  cm  anharmonicity,  the  CO 
distribution  due  to  hot  electrons  is  essentially  identical  to 
that  predicted  without  N2-  However,  on  longer  time  scales  a 
comparison  of  Figure  15  with  N2  to  Figure  8  without  N2  shows  the 
greatly  reduced  inversion  with  added  N2  and  clearly  reflects 
the  flow  of  energy  into  low  levels  of  CO. 
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Figure  14:  Evolution  of  the  CO  vibrational  distribution  due 
to  hot  electrons  in  a  pulsed  CO/N2  laser. 
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Figure  15:  Evolution  of  the  CO  vibrational  distribution  due  to 
W  processes  following  a  0.6  Msec  pulse  for  a  CO/N'^ 
laser . 


Discussion 


In  the  main  body  of  this  report  we  have  described  a  « 

full  numerical  model  for  the  excitation,  relaxation  and 

spontaneous  emission  of  a  binary  molecular  gas  mixture.  * 

For  reasonable  assumptions  on  the  system  parameters  we  have 

shown  agreement  with  the  delay  times  to  maximum  gain  in  a 

pulsed  CO  laser  system.  Further  we  have  investigated  the 

sensitivity  of  these  predictions  to  the  several  assumed 

experimental  varibles.  Clearly  before  extrapolations  of  these 

predictions  can  be  made  to  other  systems,  more  experimental 

variables  need  to  be  measured  to  permit  less  flexibility  in 

the  assumptions  of  the  numberical  model.  The  measured  CO 

W  rate  coefficients  at  room  temperature  were  important  in 

the  model  calculations,  however,  the  temperature  dependence 

of  these  rates  is  still  uncertain  even  though  it  is  likely 

that  near  resonant  exchange  is  dominated  by  long  range  forces 

which  non-resonant  processes  are  controlled  by  short  range 

processes.  Also,  the  magnitude  of  low  temperature  W  exchange 

processes  needs  to  be  determined  for  other  molecular  systems, 

i.e.  N2-CO  and  N2~N2.  Although  VT  processes  were  not  important 

on  the  short  time  scales  of  interest  in  these  pulsed  calculations, 

the  scaling  of  VT  rate  coefficients  shows  great  variability 

depending  on  the  model  adopted  as  has  been  pointed  out  before  . 

The  modeling  of  CW  laser  systems  will  not  be  reliably  performed 
until  the  questions  of  scaling  with  vibrational  level  are 
resolved. 

The  general  features  of  the  electron  coupling  rates  to 
CO  and  N2  are  clearly  understood  by  reference  to  the  unstable 
negative  ion  mechanism  responsible  for  the  large  electron  -  * 

vibrational  cross  sections.  Thus,  the  electron  couples 

efficiently  to  vibrational  levels  having  energy  nearly  equal  • 

to  the  electron  affinity  of  the  unstable  ion  but  rather 
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unefficiently  to  levels  removed  from  this  energy.  This 
feature  of  the  electron-CO  coupling  is  responsible  for  the 
partial  inversion  on  levels  around  v  =  5  to  v  =  7  and  results 
in  these  transitions  reaching  maximum  gain  at  the  earliest 
time  delays. 

One  point  indicated  in  Jeffer's  and  Wiswall's  experiments 

was  the  possibility  of  an  increased  translational  and  rotational 

temperatures  due  to  the  hot  electron  pulse  which  would  decay 

to  some  value  above  ambient  on  very  short  time  scales. 

Although  the  ejqjerimental  results  appeared  to  show  an  increase 

in  the  rotational  temperature  of  the  lasing  transitions  to 

about  750 °K  with  a  decay  time  back  to  ambient  after  the 

electron  pulse  of  2  y,sec,  the  source  of  this  increase  in 

temperature  is  uncertain.  Preliminary  conclusions  suggest 

that  direct  momentum  transfer  from  the  electrons  to  the 

neutral  gas  molecules  is  insufficient  to  produce  a  significant 

temperature  increase.  Also,  the  cross  sections  for  rotational 

(27) 

excitation  are  too  small'  '  to  play  a  role  on  the  short  time 
scales  appropriate  to  the  pulsed  operation.  The  only  mechanism 
which  appears  important  is  vibrational  excitation  accompanied 
by  simultaneous  rotational  excitation.  On  the  short  time 
scales  appropriate  to  pulsed  behavior,  VT  processes  are 
negligible,  however,  the  simultaneous  rotational  excitation 
can  be  easily  degraded  into  thermal  energy.  Since  spatial 
gradients  in  the  electrons  are  expected  in^  Jeffer's  and  Wiswall's 
experiments,  the  heated  gas  can  expand  into  the  unheated  regions 
on  time  scales  characteristic  of  sound  propagation,  i.e.  ^sec, 
rather  than  being  conducted  to  the  wal'l  which  would  require 
considerably  longer  time  scales  even^under  high  pressure 
operation.  These  points  need  to  be  further  explored  in  pulsed 
CO  laser  systems  although  it  should  be  indicated  that  we  were 
unable  to  obtain  a  "reasonable"  fit  to  the  experimental  delay 
times  upon  including  a  time  varying  kinetrc  and  rotational 
temperature  in  the  code. 
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The  addition  of  N2  to  a  pulsed  CO  laser  was  seen  to 
reduce  the  predicted  gains  even  though  more  electron  energy 
was  coupled  into  the  laser  system.  This  fact  was  interpreted 
as  being  due  to  the  anharmonic  nature  of  W  processes.  Under 
current  investigation  is  the  behavior  of  pulsed  CO  laser 
systems  in  which  other  diatomics  are  added  having  a  smaller 
vibrational  spacing  than  CO.  By  the  same  reasoning  as  used 
in  the  N2-CO  case,  a  molecule  with  a  smaller  vibrational 
spacing  than  CO  would  cause  a  preferential  drain  of  vibrational 
energy  from  vibrational  levels  of  CO  most  closely  resonant 
with  the  lowest  levels  of  the  added  gas.  This  process  could 
tend  to  increase  the  inversion  and  subsequently  the  gains  in 
pulsed  operations. 

Lastly,  it  should  be  pointed  out  that  the  electron 

disturbances  present  in  molecular  laser  systems  are  not 

Maxwellian  in  shape  but  are  non-Maxwellian  in  which  the  energy 

region  of  maximum  coupling  to  the  diatomic  gases  is  under 
122) 

populated'  .  We  have  not  considered  this  effect  in  our 
calculations  due  to  the  fact  that  the  calculated  electron- 
molecule  excitation  rate  coefficients  are  essentially  independent 
of  electron  temperature  above  about  10,000°K.  Thus,  to  first 
order,  we  can  assume  a  constant  electron-molecule  rate  coefficient 
calculated  on  the  basis  of  a  Maxwellian  electron  distribution. 

When  this  is  done,  the  actual  electron  density  assumed  may  not 
correspond  to  that  in  the  experiment  due  to  the  under  population 
in  the  important  electron  energy  range.  This  effect  is  likely 
to  produce  only  a  factor  of  2  or  3  difference  in  the  assumed 
electron  density  in  the  worst  case  but  should  not  affect  any  of 
the  qualitative  features  indicated  by  this  study. 
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